Astronomy & Astrophysics manuscript no. aa 


©ESO 2011 


October 13,2011 





Solar-like oscillations in the G9.5 subgiant /3 Aquilae^ 

E. Corsaro^'^'**, F. Grundahl^, S. Leccia^, A. Bonanno^, H. Kjeldsen^ and L. Paterno^'^ 

'Department of Physics and Astronomy, Astrophysics Section, University of Catania, Via S. Sofia 78, 1-95123 Catania, Italy 
^I.N.A.F. - Astrophysical Observatory of Catania, Via S. Sofia 78, 1-95123 Catania, Italy 
'Department of Physics and Astronomy, Aarhus University, DK-8000 Aarhus C, Denmark 
''I.N.A.F. - Astronomical Observatory of Capodimonte, Salita Moiariello 16, 1-80131 Napoli, Italy 

Received ; Accepted 4 Oct 201 1 

ABSTRACT 

Context. An interesting asteroseismic target is the G9.5 IV solar-like star/3 Aql. This is an ideal target for asteroseismic investigations, 
because precise astrometric measurements are available from Hipparcos that greatly help in constraining the theoretical interpretation 
of the results. The star was observed during six nights in August 2009 by means of the high-resolution echelle spectrograph SARG 
operating with the TNG 3.58 m Italian telescope (Telescopic Nazionale Galileo) on the Canary Islands, exploiting the iodine cell 
technique. 

Aims. We present the result and the detailed analysis of high-precision radial velocity measurements, where the possibility of detect- 
ing individual p-mode frequencies for the first time and deriving their corresponding asymptotic values will be discussed. 
Methods. The Fourier analysis technique based on radial velocity time series and the fitting of asymptotic relation to the power spec- 
trum were the main tools used for the detection of the asteroseismic parameters of the star. 

Results. The time-series analysis carried out from ~ 800 collected spectra shows the typical p-mode frequency pattern with a max- 
imum centered at 416 pHz. In the frequency range 300 - 600 yuHz we identified for the first time six high signal-to-noise ratio (S/N 
> 3.5) modes with ^ = 0, 2 and 11 < n < 16 and three possible candidates for mixed modes (C = 1), although the p-mode identification 
for this type of star appears to be quite difficult owing to a substantial presence of avoided crossings. The large frequency separation 
and the surface term from the set of identified modes by means of the asymptotic relation were derived for the first time. Their values 
are Ay = 29.56 ± 0.10 pHz and e = 1.29 ± 0.04, consistent with expectations. The most likely value for the small separation is 
6vo2 = 2.55 ± 0.71 yuHz. 
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1. Introduction 



The search for solar-like oscillations (see Chnstensen-Dalsgaard 
|2004 Elsworth & Thompson 2004 for a summary) in main- 



sequence and subgiant stars showed a tremendous growth in 
the last decade (for reviews see e.g. [Bedding & Kjeldsen|2003 



|2006) , especially by means of the photometric space-based mis- 
sions CoRoT (Baglin et al. 2006 ; Mic hel et al.|2008| l and Kepler 
( [Borucki et al.|2010j|Koch et al.|2010j l. The latter in paiticular is 
presently providing an enormous amount of high-quality astero- 
seismic data (e.g. see Gilliland et al. 2010 for more details on the 
KAI - Kepler Asteroseismic Investigation). Photometric studies 
of a large number of solar-type stars are fundamental for statis- 
tical investigations of intrinsic stellar properties and for testing 
theories of stellar evolution (e.g. Chaplin et al. 201 1\ . However, 
as is known from the theory of solar-like oscillations, high- 
precision Doppler shift measurements are more effective for de- 
tecting p modes of higher angular degrees. At the present time 
the echelle spectrometers such as CORALIE, HARPS, UCLES, 
UVES, SOPHIE and SARG , attainin g high-precision radial ve- 
locity (RV) measurements (Marcy & Butler 1992), offer a way 
to detect solar-like oscillations in bright asteroseismic targets 
(e.g. see [Bedding & Kjeldsen||2008] for a review on ground- 
based campaigns). The spectroscopic approach to the detection 
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of solar-like oscillations will also be used in the ground-based 
SONG project ( |Grundahl et al.|2009[ ) in the near future. 



In this paper we report the detection of excess of power and 
perform a detailed analysis based on the radial velocity mea- 
surements for the evolved subgiant star /3 Aql (HR 7602, HD 
188512, HIP 98036). This star has V ^ 3"' 699 ± 0'".016 (UBV 
photometric measurements from |Oja[[1986[ ), spectral type G9.5 
IV ( (Gray et al.|2006| ), distance 13.70 + 0.0 4 pc deriv ed from the 
Hipparcos parallax tt = 73.00 + 0.20 mas ( [van Leeuw en 2007}, 
reff ^ 5160 + 100 K ([Luck & Heiter|2005) , log^ = 3.79 + 0.06 
dValenti & Fischer[2005| l and [Fe/H]^ -0.17 + 0.07 ( |Fuhrmann| 
20041. The radius R - 3.05 + 0.13 .Rq was measured by means 
of Long Baseline Interferometry ( [Nordgren et al.| 1999 1 that is 
listed in the CHARM2 catalogue for high angular resolution 
measurements ( Richichi et al.| [2005). An excess of power in 
the power spectrum (PS) of jS Aql data acquired with HARPS 
was already found by [Kjeldsen et al.[ ( [2008| l, who estimated a 
mean large separation Av = 30 pHz from stellar parameters. In 
Section |2] we describe the observations, the data reduction and 
the RV extraction developed for this star. The Fourier analysis 
and the mode characterization that led us to derive the astero- 
seismic parameters are presented in Section |3] In Section]?] we 
discuss the scaled mass, mode amplitude, and frequency of max- 
imum power of the star, where a comparison with expectations 
and a global list of stellar parameters also derived by means of 
the SEEK package ( [Quirion et al.|2010[ ). 
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2. Observation and data reduction 

The data were acquired with observations carried out during 
six nights (2009 August 5-11) by means of the high-resolution 
cross-dispersed echelle spectrograph SARG (|Gratton et al.|2001 1 
|Claudi et al. 2009 ) mounted on the 3.58 m Italian telescope TNG 
at the La Palma observing site. SARG operates in both single- 
object and long-slit (up to 26") observing modes and covers a 
spectral wavelength range from 370 nm up to about 1000 nm, 
with a resolution ranging from R = 29000 up to 7? - 164000. 
Our spectra were obtained at R = 164000 in the wavelength 
range 462-792 nm. The calibration iodine cell works only in 
the blue part of the spectrum (500-620 nm), which was used 
for measuring Doppler shifts. During the observing period we 
collected spectra with a signal-to-noise ratio (S/N) varying from 
150 to 300, a typical exposure of ~ 150 s and a sampling time of 
~ 190 s. A total of 828 spectra were collected with the following 
distribution over the six nights: 91, 133, 110, 146, 163, 185. The 
spectra were then reduced and calibrated in wavelength with a 
Th-Ar lamp, using standard tasks of the IRAF package facilities. 
No flat-fielding was applied to the spectra because of a degrada- 
tion of the S/N level. 



2.1. Radial velocity measurements 

The RV measurements were determined by means of the iSONG 

[2009) 1. 



1000 



code, developed for the SONG project (Grundahl et al. 



iSONG models the instrumental profile (PSF), stellar, and iodine 
cell spectra to measure Doppler shifts. The observed spectrum 
was fitted with a reconstructed one by using a convolution be- 
tween the oversampled stellar template, the very high-resolution 
iodine cell spectrum, and the measured spectrograph instrumen- 
tal profile. Essential to this process are the template spectra of 
j3 Aql taken with the iodine cell removed from the beam, and 
the iodine cell itself superimposed on a rapidly rotating B-type 
star, the same for all the measurements. The B-type star spectra 
allow us to construct the instrumental profile, while the /3 Aql 
spectra acquired without the iodine cell create the original stel- 
lar template when deconvolved with the PSF (see IButler et al. 



1 1996 for a detailed explanation of the method). The velocities 
were corrected to the solar system barycenter (Hrudkova et al. 
2006) and no other corrections, such as decorrelation or filtering 
by removing polynomial fits to the time-series, were applied. 

iSONG provides an estimate of the uncertainty in the veloc- 
ity measurements, cr,; these values were derived from the scatter 
of the velocities measured from many (^ 650), small (^ 1 .8 A) 
segments (chunks) of the echelle spectrum. To include the accu- 
racy of the measurements in a weighted Fourier analysis of the 
data, we firstly verified that these cr,- reflected the noise prop- 
erties of the RV measurements and their Fourier transform, fol- 



lowing a slightly modified ,Butler et al. ( 2004) approach. First we 



considered that the variance deduced from the noise level cr^, 
in the amplitude spectrum has to satisfy the relationship 



cr: 



amp 



2 



(1) 



according to |Kjeldsen & FrandseJ ( | 19921 1 and [Kjeldsen &| 
Bedding ( 1995 1. This procedm^e yields new RV uncertainties 



scaled down by a factor of 1.51. Then we performed the fol- 
lowing three steps, (i) The high-frequency noise in the power 
spectrum (PS), well beyond the stellar signal (> 1000 //Hz), re- 
flects the properties of the noise in the RV data, and because we 
expected that the oscillation signal is the dominant cause of vari- 




3.0 
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Fig. 1. Upper panel: cumulative histograms of |r,/cr,| for SARG 
data. The diamonds show the observed data, and the solid 
curve shows the result expected for a Gaussian-distributed noise. 
Lower panel: ratio / of the observed to the expected histograms, 
indicating the fraction of "good" data points. An excess of out- 
liers is evident for |r,7cr,| > 1.5. 



ations in the velocity time series, we need to remove it to analyze 
the noise. This we made iteratively by finding the strongest peak 
in the PS of the velocity time-series and subtracting the corre- 



sponding sinusoid from the time-series data (see Frandsen et al. 
1995, Section 4.3). This procedure, namely the pre-whitening, 
was carried out for the strongest peaks in the oscillation spec- 
trum in the frequency range 0-1.5 mHz, until the spectral leak- 
age into high frequencies from the remaining power was neg- 
ligible (see also Leccia et al.||2007 Bonanno et al.|2008 i. This 
left us with a time series of residual velocities, r,, that reflects 
the noise properties of the measurements, (ii) We then analyzed 
the ratio r,/cr,, expected to be Gaussian-distributed, so that the 
outliers correspond to the data points that exceed the given dis- 
tribution. The cumulative histogram of the residuals is shown 
as a solid red curve in the upper panel of Fig. [T] indicating the 
fraction of "good" data points. An excess of outliers is evident 
for |r,7cr,| > 1.5. In this step we used the theoretical Gaussian 
function for the cumulative histogram, given by the expression 



F(xi) = [1 -erf{xi/xo)]. 



where 

erf(x) 



2 



- e' dt 
n Jo 



(2) 



(3) 



is the error function, is the total number of data points, 
M - max [1 - erf {xil X{))\, and - |r,7cr,|. The zero point fixed 
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Fig. 2. Radial velocity measurements of yS Aql for the entire time of observation, obtained with the iSONG code from the S ARG 
spectra (lower panel) and their corresponding noise-scaled and outliers-corrected uncertainties cr, (upper panel). 
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Fig. 3. Detail of the fifth night of observation of radial velocity measurements obtained with the iSONG code from the SARG 
spectra (filled black circles). The solid blue line represents a 3.30 min wide smoothing to enhance the p-mode oscillations pattern. 



to xo = 0.89 allows us to adjust the fit for a noise-optimized dis- 
tribution, which means that the chosen distribution minimizes 
the noise level in the weighted PS (see [Bedding et al.|2007| for 
more details on different kinds of weight optimizations). The 
reason for this optimization relies on the possibility to improve 
the mode identification by increasing the S/N of the frequency 
peaks, (iii) The lower panel of Fig. [T] shows the ratio / of the 
values of the observed points to the corresponding ones of the 
cumulative distribution function, i.e. the fraction of data points 
that could be considered as "good" observations, namely those 
that are close to unity. The quantities w, = ^/icrf f) were adopted 
as weights in the computation of the weighted PS. 

The time-series of the whole data set is presented in Fig. |2] 
(lower panel) with the corresponding noise-scaled and outliers- 
corrected uncertainties cr, (upper panel). Fig.|3]shows the details 
of the oscillation observations during the fifth night, overlaid 
with a solid blue curve representing a smoothing of 3.30 min for 
enhancing the p-mode oscillations pattern. The final data point 
number of the full observation was reduced with respect to the 
number of observed spectra owing to a consistent improvement 
of the oscillation envelope, so that a total of 818 data points was 
taken into account for computing the time-series analysis. 



3. Time-series analysis 

The amplitude spectrum of the velocity time series was calcu- 
lated as a weighted least-squares fit of sinusoids (Frandsen et al. 
1995| [Arentoft et al.||r9 98 ; Bedd ing et al.|2 004; Kj eldsen et al. 
2005 1 with a weight assigned to each point according to its un- 
certainty estimate obtained from the radial velocity measure- 
ment, as explained in Section 2.1 The result is shown in Fig.|4] 
where a clear excess of power around 420 yuHz is visible, with 
the typical pattern for p-mode oscillations in a G9.5 subgiant 
star This feature is apparent in the power spectra of individual 
nights, and its frequency agrees with theoretical expectations, as 
we will discuss in Section |4] To determine the S/N of the peaks 
in the PS, we measured the noise level (Tamp in the amplitude 
spectrum in the range 1200 - 1500 juHz, far from the excess of 
power. By means of the new weights introduced above, it was re- 
duced from 14.4 cm s"' to the final value of 12.8 cm s"', which 
corresponds to a noise level in the PS of 0.02 m^ s"^. Since 
this is based on 818 measurements, we can deduce the veloc- 
ity precision on the corresponding time-scales using the relation 
o"RMS = o"amp '^N/n, as derived by Kjeldsen & Bedding ( 1995 i, 
which gives a scatter per measurement of 2.1 m s~'. 

However, particular care has to be taken with the noise eval- 
uation within the region of solar-like oscillations. Indeed, by 
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200 400 600 800 

Frequency (/U,Hz) 

Fig. 4. Power spectrum of the weighted radial velocity measurements of /? Aql extracted with the iSONG code from the SARG 
spectra. An excess of power is clearly visible, with a maximum centered at 422 //Hz. The inset shows the normalized power 
spectrum of the window function for a sine-wave signal of amplitude 1 m s sampled in the same way as the observations. 



evaluating the noise in the amplitude spectrum in the intervals 
100 - 300 fiHz and 600 - 800 /uHz, just below and above the ex- 
cess of power, we obtained the two noise levels cr^mp"^"" - 26.2 
cm s"' and cr^'^p^^ - 15.6 cm s~\ which appear to be quite 
different. We then adopted a noise decaying accordingly to a lin- 
ear trend law, ranging between the two values within the region 
200 - 700 yuHz. 



3.1. Search for a comb-like pattern 

The mode frequencies for low-degree, high radial order p-mode 
oscillations in Sun-like stars are reasonably well approximated 
by the asymptotic relation ( Tassoul|1980| l 



v(n,i 



= Avln + - + e\ - . 



l)^V02/6, 



(4) 



where n and d are integers that define the radial order and angu- 
lar degree of the mode, respectively, Av - {vn/ - Vn_i/) is the 
so-called mean large frequency separation and reflects the aver- 
age stellar density, 6vi)2 = (vn,o - Vn-1,2) is the small frequency 
separation, a quantity sensitive to the sound speed gradient near 
the core, and e is a quantity on the order of unity sensitive to 
the stratification of the surface layers. On attempting to find the 
peaks in our power spectrum that match the asymptotic relation, 
we were severely hampered by the single-site window function, 
whose power is visible in normalized units in the inset of Fig. [4 
giving an effective observation time of ~ 1.80 days. As is wel 
known, daily gaps in a time series produce aliases in the power 
spectrum at spacings +11.57 //Hz and multiples, which are dif- 
ficult to disentangle from the genuine peaks. Various methods to 
search for regular series of peaks have been discussed in the lit- 
erature, such as autocorrelation, comb response and histograms 
of frequencies. To find a starting value for the Av investigation 
we used the comb-response function method, where a comb- 
response function CR (Av) is calculated for all sensible values of 
Av (see Kjeldsen et al.|1995j for details), representing a general- 
ization of the PS of a PS and consequently allowing us to search 
for any regularity in the spectral pattern. In particular we used 
the generalized comb-response function discussed in Bonanno 
|etaL] ( |2008l l 



C/?(Av) = f] 



2n-l 



Avj PS (vo : 



: «Av) 



(5) 



so that a peak in the CR at a particular value of Av indicates 
the presence of a regular series of peaks in the PS, centered at 
i^max and having a spacing of Av/2. It differs from a correlation 
function because the product of individuals terms is considered 
rather than the sum. For actual calculations we used N - 2 but 
we checked that our result was stable for N > 2 as well. 

To reduce the uncertainties caused by noise, only peaks 
above 300 fiHz and with amplitude > 0.5 m s ' in the ampli- 
tude spectrum, corresponding to a S/N > 3.5, were used to com- 
pute the CR. We determined the local maxima of the response 
function CR (Av) in the range 8 < Av < 50 //Hz as first step. 
In this case the result was seriously affected by two very strong 
peaks at 11.57 //Hz and 23.14 //Hz, corresponding to once and 
twice the value of the daily gap, respectively, which consistently 
reduced the strength of the peak corresponding to Av. We then 
decided to restrict the search range to 26 < Av < 50 //Hz as sec- 
ond step to completely exclude the daily alias peaks. The result- 
ing cumulative comb-response function for this range, obtained 
by summing the contributions of all the response functions, had 
the most prominent peak centered at 28.90 + 0.45 //Hz as shown 
in Fig. |5] where the uncertainty was computed by considering 
the FWHM of the Gaussian used to fit the peak, and represented 
by the blue dot-dashed curve. The peak corresponding to three 
times the daily spacing is also visible in the right side of the plot, 
centered at 34.71 //Hz. 

3.2. Oscillation frequencies 

The comb-response function provided a guess of the large sep- 
aration, which was then used as a starting point to investigate 
the most reliable value that could represent the observed data. 
The investigation involved a parallel checking and trade-off be- 
tween two different methods, i.e. the folded PS and the echelle 
diagram. Firstly we computed the folded PS, namely the PS col- 
lapsed at the large separation value, for different values of Av 
subsequent to the guess number. The final result for the folding 
was computed for Av - 29.56 //Hz and is shown in Fig.|6] where 
the peaks corresponding to { = and £ = 2 ridges are marked 
by a dashed and a dotted line, respectively. The daily side-lobes 
for ±1 1.57 //Hz are also clearly visible and are marked with the 
same line-styles. The overlaid ridges represent the result of a 
least-squares fit to the asymptotic relation Q. It is noticeable 
how the £ = I ridge, marked by a dot-dashed line, does not 
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Fig. 5. Cumulative comb-response obtained as the sum of the 
individual comb-responses for each central frequency vq with 
amplitude > 0.5 m s"' (S/N > 3.5) in the amplitude spectrum. 
The maximum peak is centered at Av - 28.90 + 0.45 //Hz where 
the uncertainty is evaluated as the FWHM of the Gaussian used 
to fit the peak (blue dot-dashed curve). The second marked peak 
on the right, centered at 34.71 fiHz, represents three times the 
daily spacing. 
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Fig. 7. Echelle diagram overlaid on the amplitude spectrum with 
a colored background scale. The filled symbols (white and or- 
ange) represent the identified modes for £ = (circles), { = 1 
(triangles) and i - 2 (squares). The orange symbols are the 
frequencies shifted for the daily gap of +11.57 fiHz while the 
open symbols correspond to the original unshifted values. The 
ridges derived from the fit to the asymptotic relation Q are also 
marked. 




10 15 20 

Frequency modulo !^9.56 (^Hz) 

Fig. 6. Folded PS in normalized units in the case of Av - 29.56 
fiHz. The overlaid ridges represent the result for £ = 0,1,2 
mode degrees as derived by a least-squares fit to the asymptotic 
relation, which are marked as dotted, dashed and dot-dashed 
lines, respectively. The daily side-lobes for ±11.57 in the case 
of ^ = 0,2 are also clearly visible and marked with the same 
line-style for each mode degree. 



correspond to any peak in the folded PS. This could be caused 
by avoided crossings, although it could also be explained by a 
wrong identification of the modes, as we will discuss in more 
detail in Section [4] We note that the £ - 2 ridge appears to be 
slightly shifted with respect to the position of the maximum of 
the relative peak, a result that can be explained by the limit of 
our frequency resolution. 

The second method adopted for the investigation of the large 
separation is represented by the well known echelle diagram, 
an essential tool for frequency identification in asteroseismic 



data. The echelle diagram was computed for a list of 9 high S/N 
(> 3.5) frequencies directly obtained from the PS by means of 
the CLEAN algorithm, although a further consideration regard- 
ing this aspect is required to explain the way the final list was 
attained. In fact, two different ways of CLEANing the PS were 
adopted for this work. In the first case the algorithm was applied 
to the weighted PS for a sufficiently large number of frequencies 
(^ 50 to find all the peaks with S/N > 3) in the range 0.01 - 
700 fiHz, then the resulting values were restricted to the interval 
200 - 700 fiHz, obtaining a list of 20 frequencies. For the sec- 
ond case the PS region below 200 fiHz was at first completely 
pre-whitened. Then 20 frequencies were identified on the new 
resulting weighted PS in the range 200 - 700 fiHz. Comparing 
both lists of frequencies, the most unstable ones, i.e. those that 
were not present in both lists, were rejected and only the first 
nine frequencies were considered. These frequencies, showing 
an amplitude > 0.5 m s"' (or S/N > 3.5), were almost the same 
in both cases, however, the first frequency list was selected be- 
cause the chance of the identification of one more £ = 2 mode 
was possible for this case. In Fig. |7 we show the final result 
for the echelle diagram superimposed on a colored scale back- 
ground representing the amplitude spectrum, where the plot was 
computed in the same way as in Mathur et al. (201 Ij). The filled 



symbols (white and orange) represent the 10 identified modes 
for £ - (circles), £ = 1 (triangles) and £ - 2 (squares); the or- 
ange symbols are the corrected frequencies, i.e. those shifted for 
the daily gap of +1 1.57 //Hz from their original values (the cor- 
responding open symbols) and reported in Table [T] The aliasing 
considerably aff'ects the amplitude spectrum with the presence 
of several fictitious peaks, which appear as strong spots in the 
diagram without frequency symbols overlaid. In particular, the 
two daily side-lobes arising from the strongest mode {£ = 0) 
are clearly visible as two red spots on the right-hand part of the 
diagram. The uncertainties on frequencies are listed in Table [T] 



5 



Corsaro et al.: Solar-like oscillations in the G9.5 subgiant /3 Aquilae 



Table 1. Mode identification for p Aql, in the frequency range 
300-600 //Hz. The modes listed show an amplitude > 0.5 m 
s"' (or S/N > 3.5). The corrected frequencies reported in the 
fourth column include the frequencies shifted for the daily gap 
of +11.57 /(Hz. The fifth column represents the uncertainties as 
derived by means of the analytical relation (|6|. The ( - \ fre- 
quencies reported without any radial order are possible avoided 
crossings. 



0, n = 12 13 



i 


n 


S/N 


Corrected frequency" (//Hz) 


Uncertainty'' (//Hz) 





12 


6.4 


393.02 (+11.57) 


0.22 





13 


11.3 


422.43 


0.22 





16 


3.7 


511.20 (-11.57) 


0.34 


1 




4.0 


429.27 


0.16 


1 




4.7 


471.04 


0.20 


1 




3.5 


546.71 


0.38 


2 


11 


3.5 


389.87 


0.22 


2 


12 


3.5 


420.67 


0.21 


2 


16 


5.1 


538.03 (-11.57) 


0.42 



" The raw frequencies can be evaluated by adding the daily frequency 
reported in parentheses. 

The actual uncertainties can be several times larger up to the limit 
of the formal resolution of 2.2 //Hz, because the modes are not co- 
herent. 



Table 2. Asymptotic parameters as derived by a linear weighted 
least-squares fit to the asymptotic relation Q. 



Av (/iHz) 6vm (//Hz) 



29.56 + 0.10 2.55 + 0.71" 1.29 + 0.04 

The most likely value for the small separation was derived by us- 
ing the definition from the asymptotic relation with the couples of 
frequencies (viz.o, vu,2), (vi3,o, V12.2) and (vie o, Vi6,2 - Av). 



and were evaluated by using the analytical relation provided by 
[Montgomery & 0'Donoghue| ( [T999l ) 



cr(v) 



1 <cr,> 016 



//Hz, 



(6) 



where = 818 is the total number of data points, T - 5.21 
is the total duration of the run in days, {cr^) - 2.72 m s"' is 
the average uncertainty on radial velocity each data point and 
A is the amplitude per mode as derived in Section |3.3| This re- 
lation holds exactly for coherent oscillations, hence we remark 
that the estimated uncertainties represent only a lower limit to 
the real uncertainty value because these oscillations are not fully 
coherent. An upper limit to these uncertainties can be fixed to the 
formal frequency resolution, given as the reciprocal of the total 
duration of the run, which is 2.2 //Hz for this data set. The com- 
plete identification of the p modes is reported in Table[T]together 
with their S/N, where only three frequencies out of nine were 
shifted by the daily alias. The { - \ modes reported without any 
radial order number are potential candidates for mixed modes. 

By means of a linear weighted least-squares fit to the asymp- 
totic relation of the I - Q frequencies, the final value of 
Av - 29.56 + 0.10 //Hz was obtained, together with the con- 
stant € - 1 .29 + 0.04. The most likely value for the small sep- 
aration was derived by using the definition from the asymptotic 
relation with the frequency pairs (vi2_o, vn 2), (^13,0,^12,2) and 
(vi6,o, 1^15,2), where vi5,2 - 508.47 //Hz was computed directly 
from vi6,2 = 538.03 //Hz by adopting our value of the large 
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Fig. 8. Smoothed amplitude spectrum showing the amplitude per 
radial mode computed in the range 200 - 700 //Hz. The maxi- 
mum amplitude A^ax = 76+13 cm s"' occurs at v^^j^ - 416 
/iHz. The positions of the identified t - Q frequencies as derived 
from asymptotic relationship Q are also marked. 

separation. This led to 5vq2 - 2.55 + 0.71 //Hz, as reported in 
Table |2] but because it is comparable to the frequency resolu- 
tion, its uncertainty is relatively high (Fig.[7]i. As a consequence, 
the reliability of this result requires additional investigations and 
a longer data set. 



3.3. Mode amplitudes 

The evaluation of mode amplitudes is interesting for the discus- 
sion on the L/M scaling by extrapolating from the Sun, as we 
will discuss in Sec. [4] As is known from the theory of solar- 
like oscillations, the amplitudes of individual modes are affected 
by the stochastic nature of the excitation and damping (e.g. see 
Kjeldsen & Beddingj[1995|. To measure the oscillation ampli- 
tude per mode in a way that is independent of these effects, we 
followed the approach explained in Kjeldsen et al.| ([2005 2008J. 
This involves the following steps: (i) heavy smoothing of the 
weighted PS by convolving it with a Gaussian whose FWHM is 
fixed to 4Av (which is the value considered as a standard that al- 
lows comparisons since the amount of smoothing affects the ex- 
act height of the smoothed amplitude spectrum), which is large 
enough to produce a single hump of power that is insensitive to 
the discrete nature of the oscillation spectrum; (ii) conversion of 
the smoothed PS to power density spectrum (PSD) by multiply- 
ing by the effective length of the observing run (that is given by 
the reciprocal of the area under the spectral window in power, 
namely 6.40 //Hz for this data set); (iii) subtraction of the back- 
ground noise, which we computed as a linear trend in the interval 
200 - 700 //Hz, ranging from 26.2 cm s"' to 15.6 cm s"'; (iv) 
multiplication of the result by Av/c where c - 4.09 (which rep- 
resents the effective number of modes that fall in each segment 
of length Av as evaluated in the case of radial velocities) and 
taking the square root to convert to amplitude per oscillations 
mode. The result is shown in Fig.|8]for the range 200-700 //Hz, 
where A^^x = 76+13 cm s"' centered at v^ax = 416 //Hz, which 
is assumed to be the frequency of maximum power and agrees 
with the result of Kjeldsen et al. (2008). The uncertainty on the 
amplitude is evaluated by means of the analytical estimation re- 
lation 



(7) 
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( [JVlontgomery & 0'Donoghue|1999| l, where, as for the frequen- 
cies case, = 818 is the total number of data points and 
(cTv) is once more the average uncertainty on each data point. 
Again, the uncertainty in amplitude represents a lower limit to 
the real value. The amplitude distribution is evaluated for the ra- 
dial modes only, but the calculation of the amplitude in the case 
of ^ = 1 and { - 2 modes is straightforward, namely it can be 
obtained by multiplying the result for a factor of 1.35 and 1.02 
respectively, representing the relative strength given by the spa- 
tial response function (see |Kjeldsen et al. 2008 for more details). 
The result derived in this work is not far from the value obtained 



Table 3. Global list of stellar parameters for p Aql. 



by Kjeldsen et al. (2008 1. 



4. Discussion 

A complete discussion on the evolutionary state of this star 
goes beyond the scope of this work. Nevertheless, the identified 
modes provided a reliable estimate for the mean large separa- 
tion, as described in Section |3.2| which can be used to derive 
the scaled mass for this star according to the fact that Av scales 
approximately with the square root of the mean density of the 
star ( Cox & Smith|1981 1. From the scaling relation extrapolat- 
ing from the solar case (e.g. see Bedding et al.|2007 1 



' M 



R 

R^ 



Avr 



(8) 



where the radius is provided by Nordgren et al. ( 1999 1 and 
Avo - 134.9 iiHz, we obtained a mass of M = 1.36 + 0.17 



Mq, which agrees very well with the value found by Fuhrmannj 
( 2004| l. By considering a luminosity of L = 5.63 ± 0.16 Lq, as 
derived by means of visual magnitude (Oja 1986[ l, bolometric 
correction ([Flower 1996), Sun bolometric magnitude (Cox & 
Pilachowski|2000 1 and Hipparcos parallax (|van Leeuwen|2007]l, 



and the scaled value for the mass, we were able to compute 
the value for the amplitude by using the relation provided by 
Kjeldsen & Bedding] ([T995]l: 



L/Lo 



(23.4+1.4) cms- 



(9) 



which gives Aqsc = 97 ± 6 cm s^^ which fairly agrees with the 
value obtained in Section[33] Moreover, the expected frequency 
of maximum power can be evaluated from scaling the acous- 
tic cutoff frequency from the solar case. We computed this fre- 
quency by considering the relationship 



{RIR^Y ^|T,«|5111K 



3.05 ■ lO^Hz 



(10) 



( Kjeldsen & Bedding|1995 i, whose result gives v^^x - 472 + 72 
yuHz, consistent with the value obtained on the smoothed ampli- 
tude spectrum. Concerning the values for the mean large separa- 
tion and the frequency of maximum power it is outstanding that 



they fit the power law relation of Stello et al. ( 2009) quite well, 
which represents a considerable validation to the reUabiUty of 
our results. Moreover, an independent measure for the large sep- 
aration by using the Eq. ([8]), with the radius from Nordgren et al. 
( |1999 1 and the mass from Fuhrma nnj (|2004| l is also compatible 



with the value presented in Section 3.2 



Lastly, an important note regards the dipole frequencies re- 
ported in this work, which are expected to deviate strongly from 
their asymptotic values, especially for an evolved subgiant star 
like p Aql. Indeed, the presumed ( — \ modes identified here 
are possibly bumped because of avoided crossings. Although 



Stellar parameter 



Value 



Source 



M 
R 

<P> 

Z 

X 

Age 
L 

logg 

Av 

A 



1.36 ±0.17 Mo 
3.05 ± 0.13 

0.0676 ± 0.0004 g cm-^ 
0.015 ± 0.002 
0.70 ± 0.02 
lAy_^l Gyr 
5.63 ± 0.16 Lo 

-0.02 

5160 ± 100 K 
29.56 ±0.10 //Hz 
2.55 ±0.71 /iHz 
76 ± 13 cm S-' 



This work" 



Nordgren et al. 1999 
i his work 
SEEK 
SEEK 
SEEK 
This work 
SEEK 



(Luck&Heiter2005i 
I his work 
This work 
This work 



The corresponding parameter was derived directly in this work. 



the lack of a clear I - \ ridge in Fig. [6] could also be ex- 
plained by a wrong identification of the modes, we are confi- 
dent that the large separation is correct. The reason of our be- 
lief is that different tools for its investigation, such as the comb- 
response function, the folded PS, the echelle diagram and the 
fit to the asymptotic relation, show consistent results. In addi- 
tion, its compatibility with the independent estimation and the 
agreement with the Av - v^ax power law relation mentioned in 
the above paragraph ensure a robust derivation of the result pre- 
sented. Nonetheless, we are talking about a very difficult star, 
such as V Ind (Carrier et al.|2007 ), because in the region of the 
HR diagram to which the star belongs, avoided crossings con 
siderably hamper the p-mode identification (see Bedding|201 1 



for a summary on ground-based observations across the HR di- 
agram). Therefore, a theoretical confirmation is required before 
adopting the € - \ frequencies reported in Table [T] as real fre- 
quencies of oscillations for mixed modes, and more observations 
by a multi-site ground-based project such as SONG are required 
to firmly solve the mode identification. 

4.1. Stellar parameters 



The SEEK package ( [Quirion et al.|[20T0 1 is developed for the 
analysis of asteroseismic data from the Kepler mission and is 
able to estimate stellar parameters in a form that is statistically 
well-defined. It is based on a large grid of stellar models com- 
puted with the Aarhus Stellar Evolution Code (ASTEC), which 
allow us to derive additional stellar parameters giving as input 
astrophysical quantities such as T^s ( [Luck & Heiter|2005 1, log g 
( jValenti & Fischer|2005| ) and [Fe/g] lFi3irmann|2()04| irand the 
asteroseismic values derived in this work, i.e. the large and small 
separation and the frequency of maximum power. The output list 
of parameters for yS Aql is shown in Table [3] where the l-cr er- 
ror bars on the SEEK values are computed using the Bayesian 
evaluation of the posterior distributions. The mass, radius and 
luminosity computed by SEEK also agree with the values pre- 
sented in this work. 



5. Conclusions 

Our observations of (3 Aql show an evident excess of power in the 
PS region centered at 415 juHz, clearly very well separated from 
the low-frequency power, and with a position and amplitude that 
agree with expectations. Although consistently hampered by the 
single-site window, the comb analysis and the echelle diagram 
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show clear evidence for regularity in the peaks at the spacing ex- 
pected from the asymptotic theory. The complete identification 
of six high S/N modes for £ = 0,2 led to a well-constrained mean 
large separation of Av = 29.56 + 0.10 juHz, compatible with a 
scaled value from the Sun and the value obtained by the power 
law relation (j Stello et al.||2009 1, and to a most likely value for 
the small separation of 6vo2 - 2.55 + 0.71 yuHz, whose reliability 
has yet to be confirmed. The £ - 1 modes found are presum- 
ably mixed modes but a theoretical confirmation is needed be- 
fore adopting these values as real modes of oscillations for this 
star. Moreover, our results provide a valuable proof that oscilla- 
tions in an evolved subgiant solar-like star show amplitudes that 
scale as L/M by extrapolating from the Sun. 

This campaign of observations attained with SARG led to 
high-precision RV measurements by means of the iSONG code, 
which was used for the first time in this work. The time-series 
analysis of the given data set was able to provide for the first 
time global asteroseismic parameters and individual p modes to- 
gether with the evidence for mixed modes. Moreover, this re- 
sult will be extremely important to develop theoretical models 
for this star. Multi-site observation campaign with the SONG 
project is highly desirable in a near future. That would then al- 
low us to explore the solar-like oscillations for this target in a 
detailed way by providing a large number of identified modes. 
The asteroseismic and astrophysical parameters of this star will 
then be constrained properly, yielding a deeper understanding of 
solar-like oscillations in the very difficult region of the HR dia- 
gram to which /3 Aql belongs. 
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